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Abstract 

We derive a factorization theorem that describes an energetic hadron h fragmenting from a jet 
produced by a parton i, where the jet invariant mass is measured. The analysis yields a "fragment- 
ing jet function" ^/f (s, z) that depends on the jet invariant mass s, and on the energy fraction z of 
the fragmentation hadron. We show that can be computed in terms of perturbatively calcula- 
ble coefficients, Jij{s,z/x)^ integrated against standard non-perturbative fragmentation functions, 
D^{x). We also show that Y.h ! dzzg^{s, z) is given by the standard inclusive jet function Ji{s) 
which is perturbatively calculable in QCD. We use Soft-Collinear Effective Theory and for sim- 
plicity carry out our derivation for a process with a single jet, B — t- Xhiu, with invariant mass 
^ Aqqy). Our analysis yields a simple replacement rule that allows any factorization theo- 
rem depending on an inclusive jet function J, to be converted to a semi-inclusive process with a 
fragmenting hadron h. We apply this rule to derive factorization theorems for B — )• XK^ which 
is the fragmentation to a Kaon in 6 — t- 57, and for e'^e" — t- (dijets) -|- h with measured hemisphere 
dijet invariant masses. 
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I. INTRODUCTION 



Factorization theorems are crucial for applying QCD to hard scattering processes invol- 
ving energetic hadrons or identified jets. In single inclusive hadron production, an initial 
energetic parton i = {u,d, g,u, . . .} produces an energetic hadron h and accompanying 
hadrons X. Factorization theorems for these fragmentation processes have been derived at 
leading power for high-energy e~^e~ — )• Xh, 

da = ^da,®D^ , (1) 

i 

as well as lepton-nucleon deeply inelastic scattering, e~p — e~Xh, 

da = daij O O fj /p . (2) 

tions of perturbatively calculable hard scattering cross sections, da, with non-perturbative 
but universal fragmentation functions D^{z), and parton distributions fj/p{^)- The frag- 
mentation functions D^{z) encode information on how a parton i turns into the observed 
hadron h with a fraction z of the initial parton large momentum. Fragmentation functions 
are also often used for processes where a complete proof of factorization is still missing, such 
as high-energy hadron-hadron collision, H1H2 — ?■ hX. 

Another interesting class of hard scattering processes are those with identified jets. Ex- 
amples include dijet production e^e~ — )■ Xj^Xj^^Xg where Xj^ 2 are two jets of hadrons, and 
Xg denotes soft radiation between the jets. If we measure an inclusive event shape variable 
such as thrust, or hemisphere invariant masses, then the cross section for this dijet process 
has the leading order factorization theorem [2|-l5j 

da = ifajet J ® S2jet ■ (3) 

Here the J = J{s) are inclusive jet functions depending on a jet invariant mass variable s, 
S2jet is a soft function which gets convoluted with the Js as denoted by (g), and -ff2jct is a 
multiplicative hard coefficient. Another example of this type is the B — )■ Xuii'e decay in a 
region of phase space where X^ is jet-like 
factorization theorem for the decay rate is 



Aqcd ^ ''^Xu ^ "^s)- Here the leading order 

m 



dT = H Ji(S)S, (4) 

with a hard function H for the underlying b — uiue process, the same inclusive jet function 
J as in the previous example, and a "shape function" S which is the parton distribution for 
a 6-quark in the 5-meson in the heavy quark limit. 

In this paper we will analyze processes which combine the above two cases, namely both 
the fragmentation of a hard parton i into h and the measurement of a jet invariant mass. 
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Since this probes fragmentation at a more differential level, we expect it can teach us inter- 
esting things about the jet dynamics involved in producing h, and shed light on the relative 
roles of perturbative partonic short-distance effects and non-perturbative hadronization. We 
derive factorization theorems that depend on a new "fragmenting jet function " Q^{s,z). 
This function depends on s, the jet invariant mass variable, and on z, the ratio of the 
large light-like momenta of the fragmentation hadron and parton. Two interesting formulae 
involving will be derived. The first formula states that 

1 f dz 

Jii.s,l^) = j j^zg^{s,z,n) , (5) 

so that the inclusive jet-function can be decomposed into a sum of terms, , for fragmen- 
tation to a hadron h with m\ -C Tn\. This formula also leads to a replacement rule for 
factorization theorems, where we can take any process involving an inclusive jet function, 
and replace Ji — )■ to obtain the corresponding process with a fragmenting jet. 
The second formula states that to leading order in Aq^^j^/ s ^ 1 we have 

g^{s,z,f^) = J2 [ ^J^,{s/-,^^)D^{x,f^), (6) 

so that the fragmenting jet function can be expressed in terms of perturbatively calculable 
coefficients J'ij, together with the standard unpolarized fragmentation functions Dj{x,fi) 
renormalized in the MS scheme. 

To introduce the concept of and study its properties, we will specialize to a process 
with a single jet recoiling against leptons, namely B — )■ Xhiug. Using Soft-CoUinear Effective 
Theory (SCET) [7-10| we derive leading-order factorization formulae for B — )■ Xhiiy£ decay 
rates, in the region of phase space characterized by Aq^j-, ^ n^-'xh where the hadronic 

final state is jet-like, and where the energetic hadron h fragments from the jet. This b — )■ uiui 
process has the virtue of having a single jet whose invariant mass can be measured in a 
straightforward manner with available i?-factory data. Despite our focus on 5 — > XMu^ 
the results obtained can be immediately generalized to fragmentation in other processes 
where a jet invariant mass measurement is made. Two examples will be described. 

The paper is organized as follows. In Sec. [TTl we review the standard definition of the 
quark fragmentation function D^{z) and highlight features that are relevant for later parts 
of our analysis. Sec. [TTll is devoted to the process B — > Xhiue, including a discussion of 
kinematics in Sec. IIII A[ Results for relevant differential decay rates in terms of components 
of an appropriate hadronic tensor are given in Sec. IIIIBi Sec. IIVI contains the derivation 
of the SCET factorization formulae for B — t- Xhiui, and the definition of the "fragmenting 
jet function " Q^. In Sec. |V] we discuss the relations shown above in Eqs. ([5]) and ([6]). 
Conclusions, outlook, and the generalization to other processes are given in Sec. IVIIi 
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II. THE FRAGMENTATION FUNCTION D{z) 



Defining = (1,0,0,1) and = (1,0,0,-1), the hght-cone components of a generic 
four- vector are denoted hj a'^ = n ■ a and a~ = n ■ a where = = and n ■ n = 2. 
With we indicate the components of a'^ orthogonal to the plane spanned by and n'^. 
For energetic coUinear particles we will follow the convention where the large momentum is 
p~ and the small momentum is p~^. 

Let us consider a quark q with momentum k'^ fragmenting to an observed hadron h with 
momentum p^. In a frame where k± = 0, the hadron has p^ = z k~ and = {Pif^+ml) /p^. 
The standard unpolarized fragmentation function D'^^z) is defined as the integral over pj^ of 
the "probability distribution" that the parton i decays into the hadron h with momentum 



Ph 
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12| . see also 13l-ll6|. With the gauge choice n ■ A = 0, the unrenormalized quark 



fragmentation function has the following operator definition |11[ |: 

D'^iz) = -Jd'pi |^|^e^^-^"/2^Tr5^(0|^^(x+,0,a;x)|X/.)(X/.|V^(0)|0) 



(7) 



where ip is the quark field quantized on a;~ = 0, Ai'c = 3 is the number of colors, and the trace 
is taken over color and Dirac indices. In Eq. ([7]) the state \Xh) = \Xh{ph)) has a hadron 
h with momentum ph = {Phi'Ph)^ an average over polarizations of h is assumed. Boost 
invariance along the non-± direction implies that can only be a function of z = pj^/k~ 
and not k~ or individually. 

Performing a rotation and a boost to a frame where pfj;- = with left unchanged, k± 
becomes —p^ / and Eq. ([7]) can be written in a gauge-invariant form as [ill 



--e^^-^"/2--Tr^ {0\mx\OM\Xh){Xh\m\0)l^^^ , (8) 



where the field '^{x^) = ^(a;^,0,0_L) contains an anti path-ordered Wilson line of gluon 
fields, in a 3 representation 

•^{x^) = tjj{x^) Fexp ^ig J dsn-A'^{sn)^ . (9) 

We note that the form of Eq. (|8]) is not altered if we perform a Lorentz transformation to a 
frame where pj^ equals an arbitrary fixed reference value pj^^. In this case k± = {pj^^ —pj^)/z 
and pj^^ does not play any role due to the integrals over and x± in Eq. ([7]) . 

Our knowledge of the fragmentation functions is anchored to the use of factorization 
theorems to describe measurements of single-inclusive high-energy processes. Constraints 
are obtained by using perturbative results for the partonic hard collision as input. For 
example, writing out the complete form of Eq. ([1]) for single-inclusive e~^e~ annihilation into 
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a specific hadron h at center-of-mass energy Ecm, we have 

_^(e+e-^/iX) = V / —aiE,^,x,fi)D^iz/x,fi), (10) 

i 

where ctq is the tree level cross section for e'^e~ — )■ hadrons, ^ is the renormalization scale in 
the MS scheme, and Ci is the coefficient for the short-distance partonic process producing 
the parton i. In Eq. flTI]]) the sum includes the contributions from the different parton types, 
i = u, u, d, g . . . and the Cj's are calculable in perturbation theory, so measurements of 
da^/dz constrain D^. 

Model parameters for fragmentation functions have been extracted by fitting to cross 
section data for single charged hadron inclusive e~^e~ annihilation, including high statistics 



measurements at CERN-LEP and SLAG 17H20|. More recently, these data have been 



combined with semi-inclusive lepton-nucleon deeply inelastic and pp cross sections from 
HERMES and RHIC experiments, respectively, to perform a global analysis of pion, kaon 
and (anti-)proton fragmentation based on the factorized expressions for the relevant cross 



sections, with partonic input at next-to-leading order in QCD perturbation theory [2l|, |22| . 



see also 23| . These analyses confirm the universal nature of the fragmentation function, and. 



for the 7T~^, constrain the fragmentation model for the dominant D^{z) with uncertainties 



at the 10% level for z > 0.5 j20|]. There is less sensitivity to gluon fragmentation functions, 
and, correspondingly, these have larger uncertainties. 

Factorization theorems like the one in Eq. ffTU]) have been proven to all orders in 
at leading order in AgcoZ-E'cm for processes in which all Lorentz invariants like E^^ = 
{Pe+ + Pe-Y are large and comparable, except for particle masses The original proofs 
are based on the study of the analytic structure of Feynman diagrams and on a power- 
counting method to find the strength of infrared singularities in massless perturbation theory. 
Factorization is possible because only a limited set of regions in the space of loop and final 
state momenta contribute to leading power, namely the so-called leading regions which 
are hard, coUinear, and soft. For processes involving fragmentation, the leading regions 
contain a jet subdiagram that describes the jet in which the hadron h is observed [l|, 
see also ll|, l24j-|26| . Accordingly, the fragmentation function that can be constrained by 
applying factorization at leading power, corresponds to Eq. ([7]) only because the sum over 
X is dominated by jet-like configurations for the \Xh) states. Therefore, it is interesting 
to explore whether more can be learned about the fragmentation process when additional 
measurements are made on the accompanying jet. 

Here we consider what amounts to the simplest additional measurement, namely that of 
the jet invariant mass = {px Rather than using classic techniques we exploit 

the powerful computational framework of SCET. 
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FIG. 1: Kinematic configuration for a pion fragmenting from a jet Xtt. 



III. FRAGMENTATION FROM AN IDENTIFIED JET IN B ^ XMi? 



Consider the weak transition h — )■ ulvi measured with inclusive decays B — )■ Xud-Vi. The 
phase space region where is jet-hke plays an important role, because the experimental 
cuts which remove 6 — )■ c background most often restrict the final state to this region. 
Experimentally there is exquisite control over this process, e.g. in a large sample of events 
the neutrino momentum has been reconstructed by determining the recoil momentum of the 
-B, and the spectrum is available for the jet invariant mass rn\^ 27H29|. There has also been 
an extensive amount of theoretical work on this process based on the factorization theorem 
shown in Eq. (jll) 0, 0, |30-38|. From our perspective the nice thing about B — ^ Xuii^i is 
that it involves only a single jet, and hence provides the simplest possible framework to 
extend the factorization analysis involving jet functions to the fragmentation process we are 
interested in, where Xu — )■ {Xh)u. Here /i is a light-hadron fragmenting from a w-quark, 
with rrih ^ itlb- Without any loss of generality, we shall refer to as a pion vr for the 
following few sections, though we will return to the general notation h at the end. 



A. Kinematics 



In the B rest frame, the inclusive process B X^iuf can be completely described by 
three variables, often taken as the hadronic invariant mass rn\^, the square of the total 
leptonic momentum (with q'^ = Pi +Pp), and the charged lepton energy Ei. In the jet- 
like region a more convenient set of variables is p%^, ^f- where is aligned with 
the — 5-axis, and hence the jet-axis is along +z with p^^ = Ex^ ± \pxu\- 

With an identified hadron in the final state, B — )■ XTriui, there are three additional 
kinematic variables corresponding e.g. to the three independent components of Ptt- The 
orientation of the spatial axes will still be chosen such that pj^'^ = 0, as shown in Fig. [H In 
this frame the perpendicular component of the total lepton momentum vanishes, = 0, 
and = tubv^ with = (ra^ + n'^)/2. The six independent kinematic variables that we 
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will use to characterize this semi-inclusive process are: pj^^, p^^, E^, p~, p+, (pi, where 0£ 
denotes the azimuthal angle of the lepton with respect to the 2;- axis and = ±p^. By 
definition we have p~^^ < Px^, E^, = {p1^ +p~) /2 and pl^ = {p^^ + ml)/p~ . In this section we 
carry out a complete analysis of the kinematics, determining the phase space limits for the 
six kinematic variables without imposing any added restrictions or assumptions. (Later in 
Section [IV] we will specialize to the case of a fragmentation pion coUinear in the n direction 
with p^ < p~.) Note that rn?x^ = Px^^Px-K so for our process the measurement of the jet- 
invariant mass is a measurement of the invariant mass of all final state hadronic particles. 
Lepton masses will be neglected throughout, but the hadron mass will be kept for all 
calculations involving kinematics. 

For three of our six variables we can treat {Xtt) as a combined state \Xu), and hence 
{Px-K^PxnJ -E^} have the same limits as in the inclusive case, and are given by^ 

PXn 2 2 

In order to determine the limits for the remaining variables, let us first consider a frame 
where pxn = 0, so that it is as if we have Xu — )■ Xtt decay in the X„ rest frame. We denote 
the quantities evaluated in this frame by a '*'. Since 

+ ml-p\ 

K = TT- , (12) 

Zmxn 



the constraint p"^ > implies 



m^<E:< ^ . 13 



Furthermore, in this X^ rest frame there are no restrictions on the azimuthal angle 0* of the 
pion with respect to the 2;-axis, nor of that of the charged lepton 0^, i.e. < < In. 

The polar angle of the pion is also unconstrained: 

< ^: < TT . (14) 

Since p+* + p~* = 2E;, from Eq.(fT3!) 

9 

77) 

pt*<mxn + ^-p-* . (15) 
mxn 



Furthermore, since > 





mz 



Pt >-^- (16) 

Pit 



Table 2 in Ref. 



34| lists these limits for B — >• for the six possible orders of integration. 
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For the maximum p*~ = p*^^ = rnx-K, these hmits force p*"*" = rn^/mx-K and the pion 
travels along the z axis. (Interchanging p*"*" -H- p*~ gives the case where the pion travels 
along —z.) For p*~ = p*+ = (m^Tr + ^'i/^Xir) /'^i we have a pion traveling purely in the 
_L-plane with maximal energy. Note that Eqs. ( |T5l) and ( 1T6|) imply other limits such as 
< Ip^ILx / as well as p+* < p^J' = mxn, and 



pt* = E:-^Ef- ml cose: (17) 

for all cos^*. This holds true for both cases p^* < p^* and p^* > p~* , which correspond, 
respectively, to < cos6'* < 1 and —1 < cos6'* < 0. 

Let us now perform a boost along the 2;-axis with velocity vx-k = vx-k to the frame 
where the i?-meson decays at rest, which requires 

VXtt = ^ = — ■ — (.iOj 

Exn Pxn + Px. 

where < vx-k < 1- Boosting Eqs. f lTSj) and f lT6|) yields the final result for the phase 
space boundaries: 

f<r.<p-.., f<p:<pdi-f\,f. (19) 

PXtt Pit \ Px-nJ PXn 

Equivalently, for the opposite order of integration, 

^<pi<pi„ ^<p-,<p-Jl-4-)^^- (20) 

PXn Pt V PXttJ PXtt 

Finally, (j)} = (j)e since the boost is along the z-axis. Hence 

< 0£ < 27r . (21) 

B. DifTerential decay rates 

In this section we derive the fully differential decay rate for B — > Xiriu employing only 
the Lorentz and discrete symmetries of QCD, without dynamical considerations. We work 
in the B rest frame, and it is convenient to start by using the six independent variables: q"^, 
Ei, Ep, p^^, p^^, p^^. For the fully differential decay rate we have 

d^T TT^ A 2 

dq^ dE^ dEu dp^^ dp J dp^^ ^ (2^ 2K(27r)3 ^^^^^^^ ~ ^ ) ' (^2) 

where d'^pT,/[2ET,{2n)'^] is the phase space for the pion, E^^ = a/p^ + ml, and 

{B\Hl\Xniu) {Xniu\Hw\B) 
= 1^2^ 1^ (27r) 5 (ps - px^ -Pe- Pp) 

X l.s. 

= 16n Gl\V^t\^L^^Wai3. (23) 
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For the S-states we use the relativistic normahzation {B{p)\B{q)) = 2Eb (27r)^(5^(p — q). 
In Eq. fl23p the effective wealc Hamiltonian is 

AC 

Hw = ^ Vub [ui^^PLh] {IYPli^i) , (24) 
v2 

where Pl = {1 — 75)/2, and factoring the leptonic and hadronic parts of the matrix element 
gives the leptonic tensor L°^, and the hadronic tensor Wai^- L°'^ is computed without 
electroweak radiative corrections, so L""^ = Til^i'-f'^Pi ^p'j'^ Pl]. The hadronic tensor in the 
B rest-frame in full QCD is 



X 

I c/^xe"^"-^ ^{B\jf{x)\X'K) (X7r|J^"(0)|5) , (25) 



2m B 

^ X 



Airm 

with the flavor changing weak current = u-^^PLh. We have = Wai3{pt^,v^,q'^) and 
we will treat this tensor to all orders in a^. It can be decomposed using Lorentz invariance, 
parity, time reversal, and hermiticity into a sum of scalar functions, so 

L^f = 2 {p'^pI + p^Ip'^, - g"^ pe-p,-t e"''"" pe^ p,^) , 

= -gafi Wi + VaVp W2 - 1 ea^^u v^q" W3 + q^qp W4 + {Vaq/B + vpqa) 
+ {vaP-Kii + W^a) - ie^pf^u p^q" - i e^Pf^u v^p^ Ws + P-n^P-np 

+ {P7ra(ll3 + PTrpQa) WiQ , (26) 

with the convention 60123 = 1- The scalar functions Wi depend on the four independent 
Lorentz invariants g^, v ■ q, v -p^ and p^ ■ q, or four equivalent variables from our desired set, 

= W.ipj,^, p-^^, pt, p-) . (27) 

To derive Eq. (!27|) recall that the leptonic variable q^ equals msv^ —p^^, and can be traded 
for p^^- Also recall that = p^^^pj^^. Since the Wi do not depend on (pi or Ei we can (if 
desired) integrate over these variables without further information about the functional form 
of the Wi. In Eq. ( 126|) the Wi=i-5 are analogs of the tensor coefficients that can appear in 
the inclusive B — > X^jfz/ decay, but here they induce a more differential decay rate because 
of the identified pion. The VFj=g_io have tensor prefactors involving and have no analog 
in the inclusive decay. 

Contracting leptonic and hadronic tensors we find 

L"^ W^p = 2q^ Wi + {AEeE, - q^) W^ + 2q^ {E^ - E,) W^ 

+ (4E^Pp ■ p^ + 4Ep pi-p^- 2E^q^) W^ + 2q^ (pi ■ p^ - pp ■ p^) W; 
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+ 4 {EiPp -p^-E^pi- p^) Ws + {Api- p^Pp ■ p^ - miq') Wg , (28) 
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where W^^^^iq have dropped out since our leptons are massless. In terms of this contraction 
the fully differential decay rate is 

dq^dEidEpdpldpldp'^ 32n^ 2E^ ' ^ ^ 

where the limits on the kinematic variables are left implicit. 

We now want to express Eq. fl2^ in terms of the coordinates from the previous section: 
{Pxn^Pxn^Ee,p-,p+,(j)i}. The relations 

g2 = _ Pxn)i^B - Px^) , Ep = mB- Ee- (p^^ + pJ^)/2 , (30) 

suffice to convert the 1^1,2,3 terms. For the remaining Wi we need expressions for pi ■ p^r, 
Pp ■ Pt, and E^,. Recall that px-n = = — (P/ + Pu) is on the +f-axis, so the lep- 
tons are back-to-back in the _L-plane which is transverse to z. We perform a rotation 
about the z-axis to bring into the y-z plane with > 0. Then in spherical coordi- 
nates Pt, = {p^,pI,p^) = \pt, \ (0, sin6',r, cos6'7r), pe = Ei {sin 9e cos (j)e, sin Oi sin (j)e, cos 9e), and 
Pp = Ep (— sin 9p cos — sin 9p sin cos 9i), so 

P£ ■ Ptt = El Et, — El IjvI (cos 9i cos 9^, + sin 9^^ sin 9i sin (pi) , 
Pu ■ Pn = Ep Et, — Ep I (cos 9p cos 9t, — sin 6*^ sin 9p sin (pi) , 

E. = lip-+Pt)- (31) 

Using these expressions the two dot products can be written in terms of the desired variables. 
First note that 

\p^\cos9^ = ^{p- -p+) , (32) 

|p;|sin6'^= (IpttP - |P7rpcos^6'^)^^^ = (p^p+-m^)^^^ 

Furthermore, since pi ■ pp = EiEp — we have Pe ■ q = —Ei\(j\ cos 9i = pi + pi ■ pp = 
Ej + EiEp - ^2/2, which, with |g | = \px^\ = [E\^ - mj^^/^ = (p^^ - pJ^)/2, implies 
E.cos^, = {2El + 2EiEp-q^)/{p^^-p^J. Hence 

P i'n^B-PxJimB-pjcJ - El {2mB-Px^-Px^) 

Ep cos 9p = -{Pxn - Pxtt) - Ee cos 9g . 

Finally it is useful to note that the equality of the magnitude of the lepton transverse 
momenta implies Ei sin 9i = Ep sin 9p. 

Together the results in Eqs. f l32|) and f l33|) allow us to express L°'^Wai3 in terms of the six 
variables {Ei^p^^^p^^cpi). The only remaining ingredient needed to transform the decay 
rate to these variables is the Jacobian, which is easily derived by noting that 

= dpt dPn d<Pn = dp+ dp- dcpi . (34) 
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For the last equality we used the fact that the pion azimuthal angle becomes equivalent to 
the lepton azimuthal angle, ^0,^ d(j)£, when we rotate the pion momentum into the y-z 
plane. Although it would be interesting to consider measurements of (pe, for our purposes 
we will integrate over 0£ G [0, 27r]. Since the ^7,8,9 prefactors are linear in pi ■ or pp ■ 
they have contributions that are either independent of 0^ or linear in sin0£, and the latter 
terms drop out. In Wq the terms linear in sin0£ do not contribute and a quadratic term 
averages to J^"" d(f) ^sm"^ {(f) g) = n. All together this gives 



{KiWi + K2 W2 + K3W3 + Kq + Kj Wj 



dPx^ dPxn dp- dp+ dEi USir^ 

+ KsWs + KgWg), (35) 



multiplied by ^[(p^tt + ^-E^ — mB)(mB — Pxn ~ which gives the limits for the Ei 

integration, and 

K2 = -{Px^ - Px^)i.i^B - Pxn - 2Ee){mB - p\^ - 2Ei), 

= {Pxn - Pxn)imB ' Px^){^B " pi.)(4^^ - 2mB + Px^ + Px^), 
Re = -2 {niB - 2Ee - Px^){mB - 2Ee - [mB(p; - pt) + pt Pxn - Pn Pxn] . 

K7 = {niB - Pxn)i^B - Px7r)(4^^ - '^^B+Pxn+PXn) [^B^P" - pt)+PtPXn " P^ PXn] , 

Ks = (p; - pt)imB - PxJi^B - p\^){2mB - 4:E£ - Px^ - Pxn)^ 

Kg = - ^ + I - pt) i^B - Pxn) i.^B - P\J - 2Ee [rriB {p^ - p^) 

PXn -PXtt L >■ 

+ pt PXn - Pn P^n] }{ml{pt -P^) + 2m B {p- Ptc^ - pt Pxn) + pt PXn " Px^^ 
+ [mB{p^ - pt) + PtPx^ - Pn PXn] } 

+ 2{mB - Pxn)i^B - PxJi^B - Pxn - 2Ef){mB - ptc^ - 2Ee){ptp~ - ml] 
- (Pxn ~ Px^i^B ~ PxJiruB - p^^), (36) 

where the limits on the hadronic variables are displayed in Eqs. (fTT]) and (fT9|) . The Ki are 
useful for considering rates where the pion is observed along with a measurement of the 
charged lepton energy. 

Integrating Eq. f l5^ over the lepton energy Ei, the W3J 8 terms drop out leaving 



dPxn dPxn dp:;; dp+ 128n^ 



[Ki Wi + K2W2 + Kq We + Kg Wg^ , (37) 
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where 

Ki = {ms - PxJirnB - pjcn)iPx^ " Pxnf^ (38) 
1 



i4 



^2 = -^{Pxn-PXn 
^6 = ^ {PXn - PXnf [t^b{P^ - Pt) + PtPx^ - P^PXn] , 

^9 = Y^iPXn - PXn? {[pti^B - Pxn) + Pni^B - Pxn)T - ^^li^B - Pxn)i^B - Pxn)} ■ 

No further integrations can be performed in Eq. (137|) without first determining the hadronic 
structure functions Wi{p^^,p]^^,p:^ ,p~). 

IV. FACTORIZATION WITH A PION FRAGMENTING FROM A JET. 

Using SCET, we derive a leading order factorization theorem for the hadronic structure 
functions Wi appearing in the differential decay rates in section IIIIBI 

We focus on the region of phase space with endpoint jet-like kinematics where p~^^ ^ Pxw^ 
and with an energetic pion produced by fragmentation with p^ <^ p~. It is assumed that 
suitable phase-space cuts are applied to subtract the b ^ c background, which phenomeno- 
logically is responsible for the importance of this kinematic endpoint region. This issue is 



explored in a separate publication [39|. With p+ < p^ the boundaries for p+ and p^ in 
Eq. ( |T9l) become: 

^<P-.<P-X.. '$^Pt< {P-^ > (l - ^) + ^1 ' ^39) 

PXn Ptt [ \ PXnJ PXn) 

or, reversing the order as in Eq. (!20|) . 

^<pt< P\. , 1 ^' I ^ ^ - 1^) + # • (40) 

PXt: yP-n) \ PXn J PXw 

For jet-like final hadronic states, the relevant power counting is Exn ~ iTib, ml 3> 
"^"x ~ "^b^QCD, and p~ ~ mt- If we decompose the momentum of the remainder of 
the colhnear jet after the emission of the pion as Px = {PxyPxjPx)^ then it scales as 

~ (AQCD,"^fe, ^ymb Aqcd) = rnb{\^, 1, A) where A ~ a/Aqcd7^ is the SCET expansion 
parameter (which can be defined as A^ = Tn\^/m\ for our process). The total hadronic 
momentum p^^-w is also colhnear and scales the same way as p^. We will start by considering 
the {Xtt) system as a combined colhnear jet, to be factored from the hard dynamics at the 
scale mft, and the soft dynamics responsible for the binding of quarks in the B meson. This 
part of the computation can be carried out in SCETi with colhnear and ultrasoft (usoft) 
degrees of freedom. The energetic pion fragments from the jet and has a colhnear scaling 

~ (AQ(;.j-)/mf,, mfe, Aqcd) with much smaller invariant mass p^ <^ f^xn- The factoriza- 
tion for this second fragmentation step can be carried out by a SCETi to SCETn matching 
computation [40]. 
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We shall now consider Eq. ( 125|) at leading order in SCET to derive factorized expressions 
for the scalar structure functions Wi in the fragmentation region. We work in a frame where 
g_L = 0, which will induce a vanishing ±- label momentum for the light quark field in the 
partonic subprocess and the Xn system. Since the Xir system is coUinear, it is convenient 
to decompose momenta as p'^ = p'^ + pff where we have label momenta p^' ~ A°, p/" ~ A, 
and residual momenta p!f ~ A^. In SCET the pion phase space integral can be written as 



The same holds for the variables p'^ and Px. Thus for all the W = Wi, we can treat p'^^, 
Px-Ki P^i p^ as discrete label momenta. At the end these variables are restored to 
continuous variables using Eq. fHTl) and the analogs for phase space integrations over the 
Xtt variables. 

Matching the heavy-to-light QCD current onto SCET operators at a scale of order m^, 
at leading order one obtains |9|] 

3 

J:(x) = e^^ ™^-^' ^ ^«/°^(^) • (42) 

j = l UJ 

Here = n^V /2 + where V and V± are the (9(A°) and 0{X) label momentum opera- 
tors [lO^. The leading order SCET current becomes 

J:/'\uj) = xn,.r'<n. . (43) 

In this expression, Xn,ui = {^nWn) S,„p^, where ^„ is the n-collinear light u-quark field. The 
collinear Wilson line is defined as llOi 



W^= ^exp(-£n-A„(x)) (44) 



perms 



with collinear gluons An- Also Hv = Y'^hy, where hy is the ultrasoft heavy quark effective 
theory field, and Y(x) = Pexp (^ig J^^dsn ■ Aus{ns + x)) is a Wilson line built out of 
ultrasoft gauge fields, which results from decoupling the ultrasoft gluons from the leading- 
order collinear Lagrangian 0]. For the leading-order Dirac structures we use the basis 



41 



42] 



n = PRY, n = PRv'', n = PR , (45) 

n ■ V 

where Pr = {1 + 75)/2. Expressions of the one-loop Wilson coefficients Cj can be found in 



Ref. [7[ and the two-loop coefficients were obtained recently by several groups in Refs. 
3- 



43 
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For our derivation of the leading-order factorization formula in SCETj we follow the steps 



in Ref . 3J| , except that we will write out the dependence on V_[_ explicitly. When the current 



in Eq. fH2|) is inserted in Eq. fl25|) we have the phase 

'dS e"*^-^' e^T'-^-i^^^-^ = Sp^^.p 5v^fl / d^x e"*''-^' , (46) 



where n ■ p = nib — and we used the fact that q± = 0. Eq. fl46|) leaves discrete 5's that 
fix the label momenta, and a d'^x integration that only involves the residual momentum 
r'^ = n^r^ /2 with r+ = rrih — g+ ~ A^. Using the normalization convention in Heavy Quark 
Effective Theory, {By{k')\B^{k)) = 2v^{2TTy6^{k — k'), the leading order expression for 
Eq. ([25]) in SCET becomes 



3 



j,j' = l U],U]' 

[n,Ty^Xn,u.',o^]ix)\X7i){X7i\ [xn,uT,,n,]{0)\B,) , (47) 



X 

X 



where Xn,cu',Oj_ = S^,P So,v^ {W^ ^„) . Here Tj>^^ = 70 r],^7o and n ■ p = p is the large 
momentum of the energetic quark producing the jet. Grouping ultrasoft and coUinear fields 
by a Fierz transformation, we have 

— — ~ — — tL ~ ~ iji 



AN, 



+ ..., (48) 

where one should keep in mind that the (-8^,1 ■ ■ ■ \By) states will surround the Hv bilinear, 
and the (0| ■ ■ ■ |X7r)(X7r| • • • |0) states split the Xn,uj held bilinear into two parts. The ellipses 
in Eq. fl48p denote Dirac and color structures that vanish either because they involve an 
octet matrix T" between the color singlet \B^) states, or by parity, or because the only 
available vector for the {By \ ■ ■ ■ \By) matrix element is f^, and = 0. The form of the 
collinear product of matrix elements is pictured in Fig. |2l and the most general allowed 
parameterization is 

^ Tr 5^ ^ (0|Xn,a.',Ox (x) \Xn) (X7r|xn,..(0) |0) = 

2 6^,^>6{xA6\x^) uj I ^e-^"^'/' QZ^k^^, ^, ptp^) , (49) 



X 



where the trace is over color and Dirac indices and Eq. fl47j) implies that u = n ■ p. The 
first (5-function in Eq. f l49|) stems from label momentum conservation and the remaining 
ones from the fact that the leading collinear Lagrangian contains only the n ■ d derivative. 
The arguments of Ql^ are constrained by RPI-III invariance j47|, which requires products 
of plus-momenta and minus-momenta, or ratios of minus- (or plus-) momenta. (For our 
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FIG. 2: Sketch of the hadronic fragmentation process for B — )■ Xniu. 

case RPI-III is equivalent to invariance under boosts along the z jet axis.) The arguments 
of are also constrained by plus-momentum conservation. The light-cone variable is 
the plus-momentum of the up-quark initiating the Xn production, and at the interaction 
vertex is related to the residual (soft) plus-momentum of the 6-quark in the 5-meson by 
= — r+, as shown in Fig. [2J^ The large label partonic momentum n ■ p also is fixed in 
terms of kinematic variables: 



n ■ p = rrih — n ■ q = mf, — rriB + n ■ pxn = n 



Px.-A + o(^) . 

\ mh J 



(50) 



Since A = (9(Aqcd), the ratio p~ /Px^ is identified to leading order with p~ /u = p^ /p~ = z, 
the fragmentation variable in Eq. ([7]). 

Since p^ p~ = m\ Q'^ depends on which is non- vanishing with our choice 

of coordinates (the pion has p^ and the ±- momentum of X is —p'^)- In general it is the 
relative ±- momentum between vr, X, and the B that can not be transformed to zero. Later 
we will integrate over |p;Jr| or equivalently and study 



(51) 



which occurs in 



^Tr Jdp+Y,^(^\x 



oJx)|X7r)(X7r|x„,..(0)|0) 



The fragmenting jet function Q^^ defined in Eq. ( H9|) describes the properties of a final 
state that is collimated in the n-direction and consists of a up-quark initiated jet from which 



^ Strictly speaking the result A:+ = ^+ — r+ also encodes the presence of Wilson lines in defining these 
momenta, which ensure gauge invariance. 
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a pion fragments. Unlike the standard unpolarized parton fragmentation function DJ{z), 
^^(s, z,p^p~) carries information about the invariant mass s of the fragmenting jet and the 
direction of the fragmenting pion through p^p~- The matrix elements in Eqs. (jUj) and fl52|) 
are similar to the coUinear matrix element defining the jet function J{s) that appears in 
B Xuiu and B — )■ Xs'j. The jet function can be written as 

-i^Tr5^(0|^Xn(a;)|X,)(X„|xn,o.,ojO)|0) =5(x+)(52(xx)u;/" dk^ e''^^^' Uujk+) . 

C -tr J 

(53) 

J„(s) depends on the product s = p]^^, which is analogous to the first argument in Q^. 

The sum in Eq. ( 13^ extends over states with invariant mass up to m^^ ~ mfoAqcD, 
which are complete in the endpoint region. Hence one can write J as the imaginary part 
(or discontinuity) of a time-ordered product: 



Ju{k^uj) 



1 



Im / d'x e''-' I (0| Txn,.,ojO) Axn(x)|0> 



which is perturbatively calculable with 

Ju{s) = 5{s) + 0{as) 



(54) 



(55) 



On the other hand, involves a pion state, and therefore contains both perturbative and 
non-perturbative parts. A simple discontinuity formula like Eq. (l5^ does not exist for Q'!^ 
with the states {|X7r)}. 

Combining Eqs. f l47|) and f H9|) the hadronic tensor at leading order becomes 



Arc 



dk 



27r 



-ik+x- /2 



^ Cj,{m^,Px^)Cj{mi,,Px^)i:i 



Pv p{0) ^ p(0) 

2 ^'^ 2 



GZik^Px., 



PXn 



ptp-) {Bu\K{x)Y{i,0)K{0)\B,) 



X 



l + O 



A 



QCD 



m 



Xt. 



where 



Pv 



i + i) 



Y{x,y) =Y{x)Y\y) , i^ = n-xn^l2. 



(56) 



(57) 



48 



49| 



The matrix element of the bilocal operator in Eq. (]56l) defines the leading-order shape 
function 



^e-^^-'^/2^5,|/i„(£)F(x,0)/i,(0)|5,) = l{B,\K6{l+-zn-D)K\B,) , (58) 



with = r"*" + A;+ for Eq. fl56|) . In the limit rrib — ?■ oo, the support of / is (— oo. A]. The shape 
function accounts for non-perturbative soft dynamics in the 5-meson. Defining projectors 
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^ such that Wi = W^y P!^" , we obtain the leading power result 

Wl'^ = ^ Vx. t ' rffc^ Ql (k^Px., Vlv-. /i) + r+, /i) , (59) 

where we show explicitly the dependence on /x, the renormalization scale. Here hi = 
hi^mbyPx-jTjPxny A^) where the dependence on pj^^ is entirely from contractions in the tensors, 
while that on m^, p^^tt' comes also from loops. In terms of the Wilson coefficients, 



^ f (0) i y{0) 

2 ^'^ 2 



Pr . (60) 



The projectors P!^'^ relevant for the differential decay rates in Eqs. fl35!) and fl37|) have the 
following structure: 

pr = A, g^'' + v^'v^ + Q q'^q" + {v^q'' + v'^q^) + + {v'^p^ + t;^) + 

+G, (p^^ + g^) + z e^'^"^ qp + hi e^'^'^^pl q^ + hz e^^'^^v^p} , (61) 

where the coefficients Ai . . .Li are functions of pj^^, p'^^, pt, P^ and that are straight- 
forward to determine by inverting the result for Wa/^ in Eq. fl26|) . 
In terms of hadronic variables, Eq. ( |59l) becomes 



= -PxJ'^^dk^g:fk+Px^,P^,ptp-,f?jf{k+ + A~p^x.,f^) 
^ Jo ^ Vx. ^ 

= -Px.r^ dk+Qilk^Px^, ^,ptp-,^^S{p+x.-k^,^^) 

^ Jo ^ Vx. ^ 

= ^PxJ'^^ dk'+g:(p-xMxn-k'^), ^,PtPn, A S{k'\f,) (62) 
^ Jo \ Px^ J 



where S{p) = /(A — p) has support for p > 0. The convolution variable k'~^ = p'^^ — k~^ 
represents the plus-momentum of the light-degrees of freedom (soft gluons, quarks, and 
antiquarks) in the i?-meson, and PxniPxTT ~ k'^) is the invariant mass of collinear particles 
in the w-quark jet including the fragmentation pion. 

Evaluating the traces in Eq. fl6Up . we derive from Eq. ([35D the following factorization 
formula for the endpoint fivefold differential decay rate: 

SEo H{mB,Pxn^ Px^^ Ei, fi) p^^ 



dPxTT dPxn dp- dp+ dEi 

X dk'^ 0:(pxAPx. - k'^), ^,ptp-,ii) S{k'\ii) , (63) 
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with To = Gl |Kf>IV(1536 7r6) and 



H{mB,Px^,PxK^ Eh /^) = 4 {niB - " 2^^) \ {^3 - Px^)(2^f - + p^^ + p+J Cl 



q2 



+ 2 



C,Cs + {mB-p^^)C2Cs + 



PXu - PXn 



(64) 



The pi^- and i?£-dependence in this expression comes solely from contraction of the lep- 
tonic and hadronic tensors. The renormalized Wilson coefficients which encode hard loop 
corrections are functions Ci = Ciirrib^p'^^, jj). 

For the decay rate in Eq. fl371) which integrates over the lepton energy we obtain: 



dPxK dPxn dP^ dp-l 



ToH{mB,Px^,PxK^fJ') Pxn 



(65) 



X 



dk'^ G: (pxMxk - k'^)^ ^.ptp-,fA Sik'\ fi) , 

^ PXn ^ 



where 



H{mB,Pxn,PXn^f^) = iPxK-Pxnf i^B - Px^) {Sttib - 2px^ - Px^) 



+ {mB - PxJiPXn - PXn) Cl C2 + {Pxk - PXr. 



\2 '-"2 



+2 {niB - p\^Y Cl C3 + {niB - p]cn)iPx. ' Pxn) C2 Cs 
Hms-p^xJCl 



(66) 



The functions H and H encode contributions from hard scales, and from the kinematic 
contraction of tensor coefficients. In the phase space region where ~ Aqcd, at leading 
order in the SCET power counting, 



H{mB,Px^,Px^,f^) = H{mb,p ,0,/i). 



(67) 



The same considerations apply to the function H{mB,PxT^,PxT,, Ei, ft) in Eq. (!63|) . Most 
often it is useful to treat the Px^ dependence from the tensor contractions exactly, without 
expanding Px^ ^ Pxn^ since at lowest order in the perturbative corrections this allows [37 1 
the endpoint jet-like factorization theorem to agree with results derived in the more inclusive 



situation where 



~ ml. 



For the purpose of comparison with phenomenology, it is appropriate to derive the ex- 
pression of the decay rate which is doubly differential in the jet invariant mass and in the 
fraction z of large momentum components. Let us first integrate Eq. (l65|l over p+. At 
leading order we can set = 0, since = 0{X^). Therefore, in the chiral limit, p^ > 
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from Eq. ( 139|) . Since the pion fragments from the jet, the maximum value of is pj^^ — k'~^. 
Hence we can write 

—————— = ToH{mB,Pxn,Pxn,f^) dk^Gl{k^Vx^, z,ii)S{j)j^^-k^,ii) (68) 

where Q'^ is defined in Eq. (!5T|) . By integrating further: 
d^T r^"^ , ^ m^^ d^T 



dp' 



(69) 



where we indicate the two changes of variable explicitly. The integration boundaries are 
derived from Eq. (fTT]) . 

In Sec JVIl we will illustrate how to extract from this doubly differential decay rate infor- 
mation about the standard parton fragmentation function D'^. 

V. PROPERTIES OF g 

A. Relations with the inclusive jet function, J{s,fj,) 

If we sum over all possible hadrons h in the Xu — )■ Xh fragmentation process, then 
the fragmenting jet function can be related to the inclusive jet function Ju{s,n) which is 
completely calculable in QCD perturbation theory. Consider the equality 

\^ f d - f d + ^'^^ _ d^T 

dP^XHdPxudpldpt-d^^J^^^ 

where the sum with h G "H^ is over all final states with an identified h hadron fragmenting 
from the M-quark jet. The differential decay rate on the right-hand side involves the hadronic 
light-cone variables in the process B — )• The sum takes p'^^ p^^ and comparin g ou r 



Eq. fl68|) with the leading-order factorization theorem for inclusive B — Xuii' (see e.g. |38|) 

d^ 

dPx^ dPx^ 



lQ7r'ToH{mB,Px^,P^JPx^ / dk' + J^{p^^{p+^ - k'+), ^M)S{k' + , ^z) , (71) 

^0 



one obtains |50|: 



J2 jdzzg^ {k+p-^, z,fi)=2 {2Txf Mk+Px^.i^) , (72) 

heHu 

for j = u, where Ju is the leading-order quark jet function. The factor z under the integral 
is explained in Ref. 50|], and is necessary to provide the correct symmetry factor for states 



with identical particles. As the notation indicates, Eq. fj72|) holds for other partons j = 
{g,d,u, . . .} as well. This relation between the fragmenting jet function Qj and the jet 
function is not surprising since the set of states {\Xh)h£-Hu} is complete. The factor 2 (27r)^ 
is related to how we normalized and incorporated the phase space for h. 
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B. Relations with the standard fragmentation function Dg{x,fi) 



In the SCET notation, Eq. ([7]) can be written in terms of the coUinear g-quark field 



^P/^Z^Tr J2^{0\K,vSo,v^Xnm\Xh){Xh\xnmO) , (73) 



X 



since \pj^ \ d\pj^\ = (p^ /2) dp'^ at a fixed value of . Here /i is the MS renormalization scale. 
According to Eq. f l52p . the integral of over its first argument can be written as 



\jdptj dx^j d\^^i:iY,iW.,P^o,r^Xn{x)]\Xh){Xh\xn{m)- (74) 



27r 

2 



If we perform an operator product expansion on the right-hand-side of this equation we 
match onto a low energy matrix element that gives the fragmentation function in Eq. (173ll . 
Thus, is given by the convolution of a perturbatively calculable Jij and the standard 
parton fragmentation function. The result includes mixing between parton types: 



1 + 0,^QCD 



(75) 



where i, j = {m, c?, (7, m, . . .}. In Ref. [51[ the concept of a quark "beam function" is discussed. 
It turns out that a quark beam function is the analog of but with parton distributions in 
place of fragmentation functions (and an incoming proton in place of an outgoing pion). The 
derivation of the factorization theorem in Eq. fl75|) can be carried out in a manner analogous 
to the matching of the gluon beam function onto a gluon parton distribution, as derived in 
Ref. 52|- For the factorization theorem for the fragmenting jet function in Eq. (175|l . the 
Wilson coefficient Jij describes the formation of a final state jet with invariant mass s within 
which the nonperturbative, long-distance fragmentation process takes place. 

At tree level Eq. (175!) is easily verified. Using a free g-quark of momentum p in place of h 
in the final state in Eq. ( 175]) (and denoting the label parts of p^ by p^, pf and the residual 
parts by p^, and defining z = p~ juj), the partonic Q is 

QT^k-^u^z) = Jdx-dxU'x^e^''-^^/'^TT[^6^^.^-6,^^.m^^^^^ 

= E /^4(2vr)^5(fc+ -p,\/p7)5(p7)52(p,^)5^,^- 6,^^. 
Pi 

= 2{27if5{k+)5iuj -p~) = 2 {2Trf 5ik^uj) 5(1 - z) . (76) 
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In the second to last step we recombined the residuals and labels into the continuous 

^uip-^^Pr) — ^('^ ~V~)- The quark fragmentation function is D'^^^^i^z) = 5(1 — z). Since 
the Wilson coefficients jTij are independent of the choice of states, the tree-level coefficient 
function can be identified as 

J'^^'^k+co, z/x, fx) = 2 (271)^ 5{k+u) 5(1 - z/x) , (77) 



which satisfies Eq. (175]) . The one-loop calculation of Jij is presented in Ref. |50[ |. 



VI. Dl{z) FROM A DOUBLY DIFFERENTIAL DECAY RATE 

As a further consequence of our factorization formulae, we explore a strategy to extract 
from measurements of suitable differential 5-decay rates the standard pion fragmentation 
function -0^(2;) for values of z that are not too small, such as ^ > 0.5. Ultimately we 
anticipate that fragmenting jet functions will be useful for many other processes (including 
hadron-hadron collisions) for which a factorization theorem like Eq.f l65p can be derived 
involving Q^. The phenomenology of i?-decays is particularly instructive to this purpose 
since it allows to concentrate on single jet production avoiding the kinematical complications 
of more involved scattering processes. 

For comparison with phenomenology, we are interested in the doubly differential decay 
rate in Eq. (l69|l . We aim at writing a factorization formula of the type: 



Eo / — Huj [mb,mj^^,-,fi] D]{x,fi) (78) 



dm\^ dz —r- 

where H^j is calculable in perturbation theory and the cut refers to a suitable interval in 
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over which we integrate. We shall argue that in the "shape function OPE" regime 
{v^Tt Pxtt ^QCd) it is possible to write a factorization formula involving D^, which 
does not specify the invariant mass of the final-state jet. 



According to the discussion in Ref. [38|, the shape function can be written as a convo- 
lution when integrated over a large enough interval [0, A] such that perturbation theory is 
applicable at the scale A: 

/>oo 

S{u)= dco'Co{co-uj')F{co') , (79) 
Jo 

where Co is the 6-quark matrix element of the shape function operator calculated in pertur- 
bation theory and F is a non-perturbative function that can be determined by comparison 
with data. F falls off exponentially for large u' and all its moments exist without a cutoff. 
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Combining Eqs. ( 168|) . (169|) and (175|) . the integration over pj^^^ leads to 



X / dk+J^Ak+^,-,fi] / duj'Co{pt^-k+-u',fi)F{uj') (80) 

^0 \ PXn ^ /Jo 



if Px'^^'^ ^ Aqcd- Let us now perform a Taylor expansion of the perturbative kernel Co 
around u' = 0: 



CoiP^. -k^- = Coipjc. - k^) - ^' C'oiP^c. -k^) + ... (81) 



Since 



38| 



poo poo 

/ duj'F{Lo') = l and / c/w' w'" F(w') = C'(Aqcd) > (§2) 
Jo Jo 

the /^^-convolution integral in Eg. (1801) can be written as 

\ dk^ Ju, k+ -, /i Co(pi, - A;+, /i) + . . . (83) 

where the dots indicate terms suppressed by increasing powers of ^qcd /Px-k phase- 
space region where the jet becomes less coUimated and increases its invariant mass {pxtt ^ 
Aqcd)- Hence, for PxT^ > Aqcd and p^™'"'' < p^^, at leading order in the Aqcd/Px-k- 
expansion, 



j=u,u,d,a- JrnlJniB VPXiv) V PXn / 

X / ciA;+ Ju, k+ -, fi Co(p+, - k+, /i) . (84) 

^0 \ PXn ^ / 



By identifying 



^ ^ ^ Jmj,Jmn [PxJ V / 

X / dk^ Ju, k+ -, /X Co(p+, - k+, fi) , (85) 

^0 V Pxn ^ J 

we see that Eq. flTHjl is satisfied at leading order in SCET. Note that to obtain the complete 
inclusive Huj there are additional hard corrections from processes beyond those treated in 
the jet-like region, so Eq. ( 185|) does not give the complete expression for Huj. Following the 
same steps, one can also test the consistency with the factorized expression 

dml^dzdE, = X ~ y^'' X ' ^) ""^^'^ ' ^^^^ 

where H{mB,Pxn = ^Xtt/Pxtt^Pxtt^ ^e, IJ-) in Eq. (JM]) replaces H in Eq. ([85]). 
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VII. CONCLUSIONS 



Using Soft-Collinear Effective Theory, we have derived leading-order factorization formu- 
lae for differential decay rates in the process B — )■ Xhiv where /i is a light, energetic hadron 
fragmenting from a measured M-quark jet. We obtained results for differential decay rates 
with various kinematic variables, for example 



ST 



fPxh 

dPxhdPxhdPh Jo 

where Fq is a constant prefactor, H encodes contributions from hard scales, z = p^/Pxh 
and S is the leading-order shape function. is the novel leading-order fragmenting jet 
function: at variance with the standard parton fragmentation function Df{z), it incorporates 
information about the invariant mass of the jet from which the detected hadron fragments. 

We have also shown that it is possible to extract D^{z) from a suitable B — )■ Xhiu 
differential decay rate, for values of z that are not too small, like z > 0.5. 

Moreover, our analysis implies that to obtain a factorization theorem for a semi-inclusive 
process where the hadron h fragments from a jet, it is sufficient to take the factorization 
theorem for the corresponding inclusive case and make the replacement 

JAk^u:) g^ik+u, z) dz , (88) 

where Jj is the inclusive jet function for parton j and the additional phase space variable 
is z = Ph/p^hi momentum fraction of the hadron relative to the total large momentum 
of the Xh system. This replacement rule is consistent with integration over the phase 
space for h. Applying Eq. fl88l) to the factorization theorem given schematically in Eq. (jl]) 
we derive the following factorization formula for the doubly differential decay rate in the 
process B — )• XKj: 

dE \z ^ ^{2Tif ^'^^^K' / dk+gf{k^mb,z,fi)S{p1^K-k+,^) (89) 
= ^^^^Hsip'^K.n) rf/c+^ Jsj{k^mb,^,^?jDf{x,n)S{p\K-k^,n) , 



with pj^^ = niB — ^E^, and the Tqs and Hg are defined in Eq. (5) and Eq. (Al) of Ref. [38 1. 
The soft function S is the same one as in endpoint B — )■ X^'j. The jet Wilson coefficients 
J'sj are process independent and calculable in perturbation theory. At tree level: 

J^l^%k+uj, z/x, fi) = 2 (27r)3 6{k+uj) 6{1 - z/x) . (90) 

Analogously, for the process e^e^ — ?■ (dijets) + /i we apply Eq. flHHj) to the factorization 
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theorem for e'^e — )■ (dijets) in Eq. (|3]) to obtain the factorized differential cross-section 



where cxo is the tree level total cross-section which acts as a normalization factor, Q is the 
center-of-mass energy, and are hemisphere invariant masses for the two hemispheres 
perpendicular to the dijet thrust axis. Since here we assume that it is not known whether 
the hadron h fragmented from the quark or antiquark initiated jet, we have a sum over 
both possibilities in the factorization theorem. For the definitions of ctq, -ff2jet, and S'2jet see 
Ref. |j] whose notation we have followed. 

The factorization formulae derived with our analysis should allow improved constraints 
on parton fragmentation functions to light hadrons, by allowing improved control over the 
fragmentation environment with the invariant mass measurement, as well as opening up 
avenues for fragmentation functions to be measured in new processes, such as -B-decays. 
We also expect that further study based on the definition of the fragmenting jet function, 
will contribute to a better understanding of the relative roles of perturbative partonic short- 
distance effects and non-perturbative hadronization in shaping jet properties and features. 

VIII. ACKNOWLEDGMENTS 

This work was supported in part by the Office of Nuclear Physics of the U.S. Department 
of Energy under the Contract DE-FG02-94ER40818, and by the Alexander von Humboldt 
foundation through a Feodor Lynen Fellowship (M.P.) and a Friedrich Wilhelm Bessel award 
(I.S.). We acknowledge discussions with F. D'Eramo, A. Jain, and W. Waalewijn. 



[1] J. C. Collins, D. E. Soper, and G. Sterman, Adv. Ser. Direct. High Energy Phys. 5, 1 (1988), 
| hep-ph /0409313|. 

[2] G. P. Korchemsky and G. Sterman, Nucl. Phys. B 555, 335 (1999), [hep-ph /9902341 1. 
[3] G. P. Korchemsky and S. Tafat, JHEP 10, 010 (2000), |hep-ph/0007005|. 
[4] S. Fleming, A. H. Hoang, S. Mantry, and I. W. Stewart, Phys. Rev. D 77, 074010 (2008), 
|hep-ph/0703207]. 

[5] M. D. Schwartz, Phys. Rev. D 77, 014026 (2008), [arXiv:0709.2709] . 



dM^ dM^ dz 





(91) 



24 



G. P. Korchemsky and G. Sterman, Phys. Lett. B340, 96 (1994). 

C. W. Bauer, D. Pirjol, and I. W. Stewart, Phys. Rev. D 65, 054022 (2002), |hep-ph/0109045]. 
C. W. Bauer, S. Fleming, and M. E. Luke, Phys. Rev. D 63, 014006 (2001), |hep-ph/0005275|. 
C. W. Bauer, S. Fleming, D. Pirjol, and L W. Stewart, Phys. Rev. D 63, 114020 (2001), 
|hep-ph/0011336|. 

C. W. Bauer and I. W. Stewart, Phys. Lett. B 516, 134 (2001), [hep-ph/010700T] . 
J. C. Collins and D. E. Soper, Nucl. Phys. B194, 445 (1982). 
J. C. Collins, Nucl. Phys. B396, 161 (1993), [hep-ph/9208213] . 

A. H. Mueller, Phys. Rev. D18, 3705 (1978). 

H. Georgi and H. D. Politzer, Nucl. Phys. B136, 445 (1978). 

R. K. Ellis, H. Georgi, M. Machacek, H. D. Politzer, and G. G. Ross, Nucl. Phys. B152, 285 
(1979). 

G. Curci, W. Furmanski, and R. Petronzio, Nucl. Phys. B175, 27 (1980). 
S. Kretzer, Phys. Rev. D62, 054001 (2000), |hep-ph/0003177|. 

B. A. Kniehl, G. Kramer, and B. Potter, Nucl. Phys. B582, 514 (2000), 
S. Albino, B. A. Kniehl, and G. Kramer, Nucl. Phys. B725, 181 (2005), 



hep-ph/0010289 



hep-ph/0502188 



M. Hirai, S. Kumano, T. H. Nagai, and K. Sudoh, Phys. Rev. D75, 094009 (2007), 
|hep-ph/0702250|. 

D. de Florian, R. Sassot, and M. Stratmann, Phys. Rev. D75, 114010 (2007), 
| hep-ph/0 703242|. 

D. de Florian, R. Sassot, and M. Stratmann, Phys. Rev. D76, 074033 (2007), 
[arXiv:0707.1506j. 

S. Albino, B. A. Kniehl, and G. Kramer, Nucl. Phys. B803, 42 (2008), [ 0803.2768| . 
A. H. Mueller, Phys. Rev. D9, 963 (1974). 

C. G. Callan, Jr. and M. L. Goldberger, Phys. Rev. Dll, 1553 (1975). 
J. C. Collins and D. E. Soper, Nucl. Phys. B193, 381 (1981). 
BABAR Collaboration, B. Aubert et al, hep-ex/0408068, 

BABAR Collaboration, K. Tackmann, Eur. Phys. J. A38, 137 (2008), [0 801.2985] . 
Belle Collaboration, P. Urquijo et ai, 0907.0379, 



A. K. Leibovich, I. Low, and I. Z. Rothstein, Phys. Lett. B486, 86 (2000), [hep-ph/0005124j. 



C. W. Bauer, M. Luke, and T. Manuel, Phys. Lett. B543, 261 (2002), [hep-ph/ 0205150| 
C. W. Bauer and A. V. Manohar, Phys. Rev. D70, 034024 (2004), |heirph/0312109|. 
S. W. Bosch, B. O. Lange, M. Neubert, and G. Paz, Nucl. Phys. B699, 335 (2004), 
|hep -ph/0402094 |. 

K. S. M. Lee and I. W. Stewart, Nucl. Phys. B721, 325 (2005), |hep-ph/0409 045|. 

M. Beneke, F. Campanario, T. Manuel, and B. D. Pecjak, JHEP 06, 071 (2005), 

|hep-ph/0411395|. 

B. O. Lange, M. Neubert, and G. Paz, Phys. Rev. D72, 073006 (2005), |' hep-ph/050407T] . 
F. J. Tackmann, Phys. Rev. D72, 034036 (2005), [hep-ph/0503095] . 

Z. Ligeti, I. W. Stewart, and F. J. Tackmann, Phys. Rev. D 78, 114014 (2008), 
|arXiv:0807. 1926j. 

M. Procura, I. W. Stewart, and F. J. Tackmann, In preparation. 

C. W. Bauer, D. Pirjol, and I. W. Stewart, Phys. Rev. D 67, 071502 (2003), |h ep-ph/02H069] . 

D. Pirjol and I. W. Stewart, Phys. Rev. D67, 094005 (2003). 

J. Chay and C. Kim, Phys. Rev. D65, 114016 (2002), [hep-ph/0201197] . 
R. Bonciani and A. Ferroglia, JHEP 11, 065 (2008), [0809.4687] . 



25 



[44] H. M. Asatrian, C. Greub, and B. D. Pecjak, Phys. Rev. D78, 114028 (2008), [0810 .0987J. 

[45] M. Beneke, T. Huber, and X. Q. Li, Nucl. Phys. B811, 77 (2009), [0810.1230J. 

[46] G. Bell, Nucl. Phys. B812, 264 (2009), [0810.5695J. 

[47] A. V. Manohar, T. Mehen, D. Pirjol, and I. W. Stewart, Phys. Lett. B539, 59 (2002). 

[48] M. Neubert, Phys. Rev. D49, 4623 (1994), [hep-ph/9312311]. 

[49] I. 1. Y. Bigi, M. A. Shifman, N. G. Uraltsev, and A. I. Vainshtein, Int. J. Mod. Phys. A9, 

2467 (1994), [hep-ph/9312359|. 

[50] A. Jain, M. Procura, and W. J. Waalewijn, arXiv: 1101 .4953. 

[51] I. W. Stewart, F. J. Tackmann, and W. J. Waalewijn, arXiv:0910.0467. 

[52] S. Fleming, A. K. Leibovich, and T. Mehen, Phys. Rev. D74, 114004 (2006), [ hep-ph/0607T2T] . 



26 



